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Abstract. We review a general class of models for self-organized dynamics based on align- 
ment. The dynamics of such systems is governed solely by interactions among individuals or 
"agents", with the tendency to adjust to their 'environmental averages'. This, in turn, leads 
to the formation of clusters, e.g., colonies of ants, flocks of birds, parties of people, etc. A 
natural question which arises in this context is to understand when and how clusters emerge 
through the self- alignment of agents, and what type of "rules of engagement" influence the 
formation of such clusters. Of particular interest to us are cases in which the self-organized 
behavior tends to concentrate into one cluster, reflecting a consensus of opinions, flocking or 
concentration of other positions intrinsic to the dynamics. 

Many standard models for self-organized dynamics in social, biological and physical science 
assume that the intensity of alignment increases as agents get closer, reflecting a common ten- 
dency to align with those who think or act alike. Moreover, "Similarity breeds connection," 
reflects our intuition that increasing the intensity of alignment as the difference of positions 
decreases, is more likely to lead to a consensus. We argue here that the converse is true: when 
the dynamics is driven by local interactions, it is more likely to approach a consensus when the 
interactions among agents increase as a function of their difference in position. Heterophily 
— the tendency to bond more with those who are different rather than with those who are 
similar, plays a decisive role in the process of clustering. We point out that the number of 
clusters in heterophilious dynamics decreases as the heterophily dependence among agents 
increases. In particular, sufficiently strong heterophilious interactions enhance consensus. 
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1. Introduction 

Nature and human societies offer many examples of self-organized behavior. Ants form 
colonies to coordinate the construction of a new nest, birds form flocks which fly in the same 
direction, and human crowd form parties to reach a consensus when choosing a leader. The 
self-organized aspect of such systems is their dynamics, governed solely by interactions among 
its individuals or "agents" , which tend to cluster into colonies, flocks, parties, etc. A natural 
question which arises in this context is to understand when and how clusters emerge through 
the self-interactions of agents, and what type of "rules of engagement" influence the formation 
of such clusters. Of particular interest to us are cases in which the self-organized behavior tends 
to concentrate into one cluster, reflecting a consensus of opinions, flocking or concentration 
around other positions intrinsic to the self-organized dynamics. Generically, we will refer to 
this process as concentration around an emerging consensus. 

Many models have been introduced to appraise the emergence of consensus [5| l2T| 123 } [28| 13?} 
[52| [63] . The starting point for our discussion is a general framework which embed several type 
of models describing self-organized dynamics. We consider the evolution of N agents, each of 
which is identified by its "position" Pi(t) E The position Pi(t) may account for opinion, 
velocity, or other attributes of agent "z" at time t. Each agent adjusts its position according 
to the position of his neighbors: 

(1.1) j^pi = a ^2a ij (p j - pi), dij > 0. 

Here, a > is a scaling parameter and the coefficients quantify the strength of interaction 
between agents i and j: the larger is, the more weight is given to agent j to align itself 
with agent z, based on the difference of their positions — pj. The underlying fundamental 
assumption here is that agents do not react to the position of others but to their differences 
relative to other agents. In particular, the a^'s themselves are allowed to depend on the 



relative differences, — pj. Indeed, we consider nonlinear models (1.1) in the sense that the 
coefficients may depend on the position, 

a ij =a ij (V(t)), V{t) := {p k (t)} k . 

This provides a rather general description for a nonlinear process of alignment. We ignore two 
other important processes involved in self-organized dynamics as advocated in the pioneering 
work of Reynolds, [55], namely, the short-range repulsion (or avoidance) and the long-range 
cohesion (or attraction), and we refer to interesting recent works driven by the balance of 
these two processes in e.g., [HI [271 EHl HS (49j [61] . Our purpose here is to shed light on the role 
of mid-range alignment which covers the important zone "trapped" between the short-range 
attraction and long-range repulsion. 

We distinguish between two main classes of self- alignment models. In the global case, the 
rules of engagement are such that every agent is influenced by every other agent, > 77 > 0. 
The dynamics in this case is driven by global interactions. We have a fairly good understanding 
of the large time dynamics of such models; an incomplete list of recent works in this direction 
includes [U [91 [231 E3 EE B2 SOI S3 [52], and the references therein. Global interactions which 
are sufficiently strong lead to unconditional consensus in the sense that all initial configurations 
of agents concentrate around an emerging limit state, the "consensus" p°°, 
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The first part of the paper, section [2j contains an overview of the concentration dynamics in 
such global models, from the perspective of the general framework of (1.1). 

In more realistic models, however, interactions between agents are limited to their local 
neighbors, [U O |22l 155] . The behavior of local models where some of the may vanish, 
requires a more intricate analysis. In the general scenario for such local models, discussed 
in section [3j agents tend to concentrate into one or more separate clusters. The particular 
case in which agents concentrate into one cluster, that is the emergence of a consensus or a 
flock, depends on the propagation of uniform connectivity of the underling (weighted) graph 
associated with the adjacency matrix, {a^}. This issue is explored in section [4] where we 
show that connectivity implies consensus. Thus, the question of consensus for local models 
is resolved here in terms of persistence of connectivity over time. However, even if the initial 
configuration is assumed connected, there is still a possibility of losing connectivity as the a^-'s 
may vary in time together with the positions V{t). This leaves open the main question of 
tracing the propagation of connectivity in time. 

Many standard models for self-organized dynamics in social, biological and physical science 
assume that the dependence of decreases as a function of |p* — Pj|, where | • | is a problem- 
dependent proper metric to measure a difference of positions, opinions, etc. The statement that 
"Birds of feature flock together" reflects a common tendency to align with those who think or 
act alike, jlQl EJ HH] . Moreover, "Similarity breeds connection," reflects the intuitive scenarios 
in which the influence coefficients increase as the difference of positions |p^ — pj| decreases: 
the more the a^-'s increase, the more likely it is to lead to a consensus. But in fact, we argue 
here that the converse is true: for a self-organized dynamics driven by local interactions, it is 
more likely to approach a consensus when the interaction among agents increases as a function 
of their difference |p^ — pj|. Heterophily — the tendency to bond more with the different 
rather than with those who are similar, plays a decisive role in the clustering of (1.1). The 
consensus in heterophilious dynamics is explored in the second part of the paper, in terms of 
local interactions of the form, — 0(|pi — Pj|), where </>(•) is a compactly supported influence 
function which is increasing over its support. In section [5] we report our extensive numerical 
simulations which confirm the counter-intuitive phenomenon, where the number of clusters 
decreases as the heterophilious dependence increases; in particular, if <j) is increasing fast enough 
then the corresponding dynamics concentrate into one cluster, that is, heterophilious dynamics 
enhances consensus. We are not unaware that this phenomenon of enhanced consensus in the 
presence of heterophilious interactions, may have intriguing consequences in different areas 
other than social networks, e.g., global bonding in atomic scales, avoiding materials' fractures 
in mesoscopic scales, or "cloud" formations in macroscopic scales. 

In the rest of the paper, we address a few important extensions of the self- alignment models 
outlined above. These extensions are still work in progress and we by no means try to be 
comprehensive. In section [6] we turn our attention to nearest neighbor dynamics. Careful 
3D observations made by the StarFlag project, [15] [161 El > showed that interactions of birds 
are driven by topological neighborhoods, involving a fixed number of nearby birds, instead 
of geometric neighborhoods involving a fixed radius of interaction. Here we prove that in the 
simplest case of two nearest neighbor dynamics, connectivity propagates in time and consensus 
follows for influence functions which are non-decreasing on their compact support. In section 
[7] we turn our attention to fully discrete models for self- alignment. The large time evolution in 
discrete time-steps, e.g., the opinion of dynamics in [Il[5jll3], may depend on the time-step At. 
Here, we show that the semi-discrete framework for global and local self-alignment outlined 
in sections [2}|5] can be extended, mutatis mutandis, to the fully-discrete case. In particular, 
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we recover a decreasing Lyapunov functional, a fully-discrete analogue of the semi-discrete 
clustering analysis in section |4.2| Finally, in section [8] we discuss the passage from the agent- 
based description to mean-field limits as the number of agents, or "particles" tends to be large 
enough. There is a growing literature on kinetic descriptions of such models, pT2| IT3| fT^ [28], l30| 
[36] and the references therein. Here we focus our attention on the hydrodynamic descriptions 
of self-organized opinion dynamics and flocking. 



1.1. Examples of opinion dynamics and flocking. Models for self-organized dynamics 



(1.1) have appeared in a large variety of different contexts, including load balancing in com- 
puter networks, evolution of languages, gossiping, algorithms for sensor networks, emergence of 
flocks, herds, schools and other biological "clustering", pedestrian dynamics, ecological models, 
peridynamic elasticity, multi-agent robots, models for opinion dynamics, economic networks 
and more; consult [SI El El HH [20l EH [25l EH EH EH [39l SQl HD HZl [50l E31 [5H EZl [60] and the 
references therein. 



To demonstrate the general framework for self- alignment dynamics (1.1), we shall work with 
two main concrete examples. The first models opinions dynamics. In these models, N agents, 
each with vector of opinions quantified by ^ G l d , interact with each other according to 
the first-order system, 

, v d ^-^ 0(|x 7 - — xJ) 
(1.2a) — = a^aij^j - x*) , a y = — . 

Here, < <j) < 1 is the scaled influence function which acts on the "difference of opinions" , 
|xi — Xj|. The metric | • | needs to be properly interpreted, adapted to the specific context of 
the problem at hand. Another model for interaction of opinions is 

(1.2b) -^x; = aJ2 a ij(*j ~ x *) > <Hj = v^T > <fe : = <Kl x j - X *D- 



The classical Krause model for opinion dynamics [131 15] is a time-discretization of (1.2b), which 
will be discussed in section [7] below. Observe that the adjacency matrix {a^} in the first model 



(1.2a) is symmetric while in the second model, (1.2b), it is not. 

Another branch of models have been proposed to describe flocking. These are second-order 
models where the observed property is the velocity of birds, p^ i-» E which are coupled 
to their location x^ E M, d . The flocking model of Cucker and Smale (C-S) has received a 
considerable attention in recent years, [23 } l24 f 136 } 19} 135]. 

/-i q \ d / x <Kl x j ~ x *l) i d 

(1.3a) — wi = a 2_^ a ijV v j ~ v i) , = J — where — x* = v*. 

3^i 

In C-S model, alignment is carried out by isotropic averaging. In [52] we advocated a more 
realistic alignment-based model for flocking, where alignment is based on the relative influence, 



similar to (1.2b), 



(1.3b) ^-Vi = a y^dijtvj - Vi) , = with fcj := <£(|xj - x^|). 

dt j^ l^k Yik 

Again, C-S model is based on a symmetric adjacency matrix, {a^}, while symmetry is lost in 



(l-3b), i.e. a^ ^ a^. 
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time (t) time (t) 



FIGURE 1.1. Evolution in time of the consensus model for two different influence 
functions <j> (Left figure: <j>(r) = X[o,i]> Rig ht figure: <j>(r) = .lXp.i/v^ +X[i/v^,i])' B ^ 
diminishing the influence of close neighbors (Right figure), we enhance the emergence 
of a consensus. Simulations are started with the same initial condition (100 agents 
uniformly distributed on [0, 10]). 



The models for opinion and nocking dynamics (1.2), and respectively, (1.3), can be written 
in the unified form 



N 



(1.4) 



^2 a ij(pj - Pi), 

3=1 



In the opinion dynamics, p \-> x; in the nocking dynamics, p^x. The degree U{ — N in the 
symmetric models, or Oi — <Kl x i — x jl) m ^ ne non-symmetric models. The local vs. global 

behavior of these models hinges on the behavior of the influence function, <j). If the support 
of (p is large enough to cover the convex hull of V(0) = {p&(0)}& 5 then global interactions will 
yield unconditional consensus or nocking. On the other hand, if <j) is locally supported, then 
the group dynamics in (1.4) depends on the connectivity of the underlying graph, {a^-}. In 
particular, if the overall connectivity is lost over time, then each connected component may 
lead to a separate cluster. Heterophilous self-organized dynamics is characterized by a locally 
supported influence function, 0, which is increasing as a function of the mutual differences, 
<pij = 0(|x^ — Xj|). The more heterophilous the dynamics is, in the sense that its influence 
function has a steeper increase over its compact support, the more it tends to concentrate in 
the sense of approaching a smaller number of clusters. In particular, heterophilous dynamics 
is more likely to lead to a consensus as demonstrated for example, in figure [TTT] (and is further 



documented in figures 5.2 and 5.7 below). Observe that the only difference between the two 
models depicted in figure 1.1 is that the influence in the immediate neighborhood (of radius 
r < l/y/2) was decreased, from (f) = lX[o,i] ( on ^ ne l e ^) m t° — O.lxpiA/2] ^[l/v^i] ( on 
the right): this was sufficient to enhance the four-party clustering on the left to turn into a 
consensus shown on the right. 
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2. Global interactions and unconditional emergence of consensus 



In this section we derive explicit conditions for global self-organized dynamics (1.1) to con- 



centrate around an emerging consensus. Our starting point is a convexity argument which is 
valid for any adjacency matrix A = {a^}, whether symmetric or not. We begin by noting 
without loss of generality, that A may be assumed to be row-stochastic, 



(2.1) 



1, 



1, 



,N. 



Indeed, by rescaling a if necessary we have Ylj^i a ij — 1> an d (2.1) holds when we set an :- 
1 — J2jyti a ij — 0- We rewrite (1.1) in the form 



N 



(2.2) 



dr 



3=1 



Thus, if we let Q(t) denote the convex hull of the properties {pk}k, then according to (2.2), is 
relaxing to the average value E Q(t), while the boundary of Q is a barrier for the dynamics. 
It follows that the positions in the general self-organized model ( |1.1[ ) remain bounded. 

Proposition 2.1. The convex hull of p(t) is decreasing in time in the sense that the convex 
hull, ft(t) := Conv ({Pi(£)}ie[i,iV]); satisfies 

(2.3) fi(t 2 ) C fi(ti), t 2 > *i > 0. 

Moreover, we have 



(2.4) 



max | pi (t) | < max | pi (0) | . 



Proof. We verify ( |2.4[ ) for a general vector norm | • | which we characterize in terms of its dual 
|w|* — sup w _^ (w, z)/|z| so that |p| — sup(p, w)/|w|*. Let w — w(t) denotes the maximal 
dual vector of Pi(t), so that (pi,w) = |pi|, then 

(Pi,w> = a((p i ,w> - (pi,w)) < aflpj - |pi|). 

Since (p^,w) < we have 



| Pi (t)| = (p i? w) + (p„w) < a(\pi(t)\ - |pi(t)|), 



and finally, |Pi(t)| < max^ |pi(t)| yields (2.4). 



□ 



Figure 2.1. The convex hull of the positions p^ 
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Remark. Since the models of opinion dynamics and flocking dynamics (1.4) are translation 
invariant in the sense of admitting the family of solutions {pi — c}, then for any fixed state c, 
proposition |2.1| implies 

max \pi(t) — c| < max \pi(0) — c|. 

i i 

Consensus and nocking are achieved when the decreasing Q(t) shrinks to a limit point 



m ^ {p°°}, 



max \pi(t) - p°°| 0. 



There are various approaches, not unrelated, to derive conditions which ensure unconditional 
consensus or flocking. We shall mention two: an L°° contraction argument and an I? energy 
method based on spectral analysis. 



2.1. An L°° approach: contraction of diameters. Proposition 2.1 tells us that {pi(t)}i 
remain uniformly bounded and the diameter, max^- \pi(t) — Pj(t)\, is non-increasing in time. 
In order to have concentration, however, we need to verify that the diameter of p(t) decays to 
zero. The next proposition quantifies this decay rate. 



Theorem 2.2. Consider the self- organized model (1.1) with a raw stochastic adjacency matrix 
A, pi] ). Let 

[p] := max|p^ - pj\ 

denote the diameter of the position vector p. Then the diameter satisfies the concentration 
estimate 

(2.5) - ~ a ^MV{t))[v(f)], := min^min{a^,a i/c }. 

/oo 
ipA(V(s)) ds — oo, 

then the agents concentrate in the sense that 

/t 
ipA(v(s)) ds — ? oo ^ lim max \pi(t) — Pj(t)\ 0. 

Moreover, if the decay of the concentration factor is slow enough in the sense that J°° exp(—a^(s))ds < 
oo, then there is an emerging consensus p°° E ^(0) ; 

POO 

(2.6b) / e~ a ^dt <oo ~» |Pi(*)-p°°| <e-^W[p(0)] for all i = 1, . . . , N. 



Remark. We note that theorem 2.2 applies to any vector norm | • |. 

Proof Set p := Ap. We begin with the following estimate, a generalization of [35] . which 
quantifies the relaxation of the p^'s towards their average values p^. 



(2.7) [Ap] < (l-^J[p], [p] =max| Pi - Pj | 



Put differently, (2.7) tells us that A is contractive in the matrix norm induced by the vector 
semi-norm [ • ]; since this bound is solely due to the convexity of the row stochastic A, we 
suppress the time-dependence of p and p. To proceed, fix any i and j which are to be chosen 
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later, and set rjk := min{a^, cijk} so that — rjk and ajk — % are non-negative. Then, for 
arbitrary w G l rf we have, 

(Pi-Pj,w) ^a ik (p k ,vr) -^2a jk {p k ,vr) 

k k 

^2(aik - Vk)(Pk, w) - ^(a^ - r]k)(p k ,w) 

k k 

< ^2(aik ~ Vk) m ax(p fc , w) - ^(a j/c - %) min(p fc , w) 



(1 - <ij) A ) ^max(p fc ,w) - rnin(p fc ,w)^ 



< (1 - ^ A )max(p /c - p^,w) < (1 - ^ A )max|p A , - p^||w|*. 

In the last step, we characterize the norm | • | by its dual |w|* = sup w _^ (w, z)/|z| so that 
(z, w) < |z| | w|*. Now, choose i and j as a maximal pair such that [p] = |p^ — p -|; we then 
have 



[P] 



(Pi-Pj,w> 
sup 1 — ^ < (1 - V A ) max |p fc - p £ \ 

w/0 



W 



and (2.7) now follows. 



Next, we rewrite (1.1) in the form, 



d 
dt 



e at p(t) = ae at Ap(t). 



Using (2.7) we obtain 
d 



d^W)] < [^ e£rf p(*)] = ^IMt)] < <*e<«{i - ^)[ P (t)], 

which proves the desired bound (|2.5[). In particular, we have 



(2.8) max|pi(t) - Pj{t)\ < exp (-a J^ A(v(s)) ds ) [p(0)] 

which proves ( 2.6a| ). Moreover, 



|Pi(*2) - Pz(*l 



< 



rt2 rt2 

/ pi{r)dr <amax / \pi(r) — pj(r)\ ds 

a/ exp (-a*(r))dr[p(0)], *(t) = / ^ (P(a) 



which tends to zero, |p^ (^2) — Pi(^i)| — > for £2 > £1 ^ 1? thanks to our assumption (2.6b). It 



follows that the limit Pi(t) — p?° exists, and hence all agents concentrate around the same 
limit position, an emerging consensus p°° E fi(0). The concentration rate estimate (2.6b) 
follows from d^8l). □ 



Theorem |2.2| relates the emergence of consensus or flocking of p = Ap — p to the behavior 
°f ^A(V(s))ds t °°5 an d to this end we seek lower-bounds on the "concentration factor" 
ipAi which are easily checkable in terms of the entries of A. This brings us to the following 
definition. 
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Definition (Active sets [52]). Fix 9 > 0. The active set, A(9), is the set of agents which 
influence every other agent "more" than 9, 

(2.9) A(0) := {j | Oij > 9 for any i}. 

Observe that since changes in time, — o,ij(V(t)), the number of agents in the active set 
A(0) is a time dependent quantity, denoted X(9) = \(0,t) := \A(9,t)\. 

The straightforward lower bound, i(j a > max^ 9 • X(9) yields 



Corollary 2.3. The diameter of the self- organized model (1.1) with a stochastic adjacency 
matrix A, (2.1) satisfies the concentration estimate 

(2.10) |[p(t)] < -a(max0 ■ \(0,t)) [p(t)]. 

In particular, the lower bound ijj A > iVmin^- aij, corresponding to 9 — min^- with A(0, t) = 
N, yields [36] 

(2.11) |p(t) - p°°| < exp ^-aN J m(s)ds^j [p(0)], m(s) := mina^(s). 
Remark. The bound (2.10) is an improvement of the "flocking" estimate J52j Lemma 3.1] 

f t [p(t)]<-a(m^e-X(e,t)f[p(t)]. 

Corollary |2.3| is a useful tool to verify consensus and flocking behavior for general adjacency 
matrices A = {a^}, whether symmetric or not. We demonstrate its application with the 
following sufficient condition for the emergence of a consensus in the opinion models (1.2). In 
either the symmetric or non-symmetric case, 

j N L«wf ^ <N 

a>ij = < > > n ? <?i = ^ - N - 

By proposition |2.1[ the diameter [x(t)] is non-increasing, yielding the lower bound 

TV 

JVo ii (P(t)) = -0(|xi(t)-x J -(t)|)> min 6(r) > minj(r), 

0"i r<[x(t)] r<[x(0)] 

which in turns implies the following exponentially fast convergence towards a consensus x°°. 



Proposition 2.4 (Unconditional consensus). Consider the models for opinion dynamics (1.2) 
with an influence function <p(r) < 1, and assume that 

(2.12) m := min <6(r) > 0. 

r<[x(0)] 

Then, there is an exponentially fast convergence towards an emerging consensus x°° ; 

(2.13) |xi(t) - x°°| < e- amt [x(0)]. 



Similar arguments apply for the flocking models ( |1.3[ ): since [v(t)] is non-increasing then 
[x(t)] < [x(0)] +t[v(0)] and hence 

TV 

JVoi j (P(t)) = -0(|x i (t)-x j (t)|)> min 6(r) > ^ min 0(r); 

Gi r<[x(t)] r<[x(0)]+t[v(0)] 
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is decreasing then we can set m(t) = 0([x(O)] + t[v(0)]) and unconditional flocking 

for sufficiently strong interaction so that f°° <p(s)ds = oo. In 



if 



follows from for corollary 



2.3 



fact, a more precise statement of flocking is summarized in the following. 



Proposition 2.5 (Unconditional flocking). Consider the flocking dynamics (1.3) with a de- 
creasing influence function <p(r) < 0(0) < 1, and assume that 



(2.14) 



/OO 
(f){s)ds 



oo; 



Then, the diameter of positions remains uniformly bounded, [x(t)] < < oo, and there is 
an exponentially fast concentration of velocities around a flocking state v°°, 



(2.15) 



< 



-amt 



[v(0)], 



m 



<t>{D 



Proof. Unlike the first-order models for consensus, the diameter in second-order flocking mod- 
els, [x(t)] 5 may increase in time. The bound stated in (2.15) places a uniform bound on 
the maximal active diameter. To derive such a bound observe that in the second order flocking 
models, the evolution of the diameter of velocities satisfies, 

|[v(t)]<-«<K[x(*)])[v(i)], 

and is coupled with the evolution of positions [x(i)]: since x = v, we have 

d 



dt 



[x(i)] < [v(t)]. 



The last two inequalities imply that the following energy functional introduced by Ha and Liu 

[x(t)] 



£{t) := [v(t)] + a / <f>{s)ds, 
Jo 



is decreasing in time, 
(2.16) 



a 



<t)] 

x(0)] 



W*<[v(0)]-[v(t)]<[v(0)]. 



This, together with our assumption (2.14) yield the existence of a finite > [x(0)] such that 
(2.17) a / 0( 5 ) ds < [v(0)] < a / <j>{s)ds. 

7[x(0)] J[x(0)] 

Thus, the active diameter of positions does not exceed [x(t)] < D^, and since (j) is assumed 
decreasing, the minimal interaction is Naij > 0([x(t)]) > (^(D^) which yields 

d -[v(t)]<-a^D OQ )[v(t)]. 



dt 1 



□ 



This concludes the proof of ( |2.15[ ). 

Remark (Global interactions). Proposition 2.4 derives an unconditional consensus under the 
assumption of global interaction, namely, according to (2.12) every agent interacts with every 
other agent as 

1 772 

a ij > ^<KI x ;- x jI) > >0- 
Similarly, the unconditional flocking stated in proposition\2.5] requires global interactions, in the 



sense of having an influence function (2.14) which is supported over the entire flock. Indeed, 
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if the influence function <j) is compactly supported, supp{(j)} — [0, R], then assumption (2.14) 
tells us that 



[v(0)] < a / (p(s) ds 
■Ax(o)] 



Hx(t)] 

but according to (2.16), a / (f)(s) ds < [v(0)] and hence the support of (j) remains larger 

J MO)] 



[x(0)] 

than the diameter of positions, R > [x(t)]. 



Proposition |2 . 5| recovers the unconditional flocking results for the C-S model obtained earlier 
using spectral analysis, t\-^ and ^-based estimates [231 E3 El E3 E3] . The derivations are 
different, yet they all required the symmetry of the C-S influence matrix, — faj/N. Here, 
we unify and generalize the results, covering both the symmetric and non-symmetric scenarios. 
In particular, we improve here the unconditional flocking result in the non-symmetric model 
obtained in [52l theorem 4.1] . Although the tools are different — notably, lack of conservation of 
momentum v i(^) m the non-symmetric case, we nevertheless end up with same condition 



(2.14) for unconditional flocking. 



2.2. Spectral analysis of symmetric models. A more precise description of the concentra- 
tion phenomena is available for models governed by symmetric influence matrices, — aji, 
such as ( |1.2a[ ) and (1.3a). Set = — (p) where (p) := l/NJ2iPi is the average (total 



momentum), which thanks to symmetry is conserved in time, (p)(t) = Eij a ij(Pi ~ Pj) — 0> 
and hence the symmetric system ( |1.1[ ) reads 

d N 

— qi(t) = a^2<iij(qj - q*), q* := Pi - (p). 

Let La := I — A denote the Laplacian matrix associated with A, with ordered eigenvalues 
= Xi(La) < \2(La) < • • • Xn(La)- The following estimate is at the heart of matter (here | • | 
denotes the usual Euclidean norm on 



(2.18) ^5^l*(*)| 2 = a J2 ai i^ = -f^^jl^J ~^| 2 

i i,j ij 

< - a \ 2 (L A )J2\^(t)\ 2 - 



The second equality is a straightforward consequence of A being symmetric; the following 
inequality follows from the Courant-Fischer characterization of the second eigenvalue of La in 
terms of vectors q orthogonal to the first eigenvector 1 = (1,1,...,1) T , 



(L A q, q) (1/2) Eij a*Mi ~ <H 



|2 



(2.19) \ 2 {L A ) = min ^ \ < 

Eq^o (q,q) Ei |q»l 2 

We end up with the following sufficient condition for the emergence of unconditional concen- 
tration. 

Theorem 2.6 (Unconditional concentration in the symmetric case). Consider the self- organized 



model (1.1),(|2.1|) with a symmetric adjacency matrix A. Then the following concentration es- 



timate holds 

(2.20) V p(t) < exp (-a J X 2 (L A{v{s)) )ds^ V p(0) , V p(t) := ^ Ip*(*) ~ 
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In particular, if the interactions remain "sufficiently strong" so that J°° \2{L^j>^)ds — oo, 
then there is convergence towards consensus Pi(t) — > p°° = (p)(0). 



To apply theorem 2.6, we need to trace effective lower bounds on \2(La)', here are two 



examples which recover our previous results in section 2.1 



Example #1 (revisiting theorem 2.2). If r is the Fiedler eigenvector associated with Ajy-i(A) 
with rll, then (|2.7l) implies 



, [i4r] [Ap] 
X N -!(A) = l — I < sup Vt 1 ^ 1 
M P ±i [p] 



We end up with the following lower bound for the Fiedler number 

X 2 (L A ) = 1 - X N -!(A) >l-(l-i/> A )> il) A . 



Thus, theorem 2.2 is recovered here as a special case of the sharp bound (2.20) in theorem 



2.6 The former has the advantage that it applies to non-symmetric models, but as remarked 
earlier, is limited to models with global interactions; the latter can address the consensus of 
local, connected models, consult [6] below. 



Example #2 (revising propositions 2.4 and 2.5). A straightforward lower bound \ 2 {La) > 
TVminaf/ recovers corollary 2.3, 



(2.21) 



V 



Pit) 



< exp 



m(s)ds ] V 



'p(O), 



m(t) 



mm < 



|x^)-x^)|), 



The characterization of concentration in theorem |2.6| is sharp in the sense that the estimate 
(2.18) is. Indeed, it is well known that positivity of the Fiedler number, \2(La) > 0, charac- 



terizes the algebraic connectivity of the graph associated with the adjacency matrix A, [32], [5Tj . 
Theorem |2 . 6| places a minimal requirement on the amount of connectivity as a necessary condi- 
tion for consensu^) There are of course many characterizations for the algebraic connectivity 
of static graphs. In the present context of self-organized dynamics (1.1), however, the dynam- 
ics of p = a(Ap — p) dictates the connectivity of A = A(V(t)), which in turn, determines the 
clustering behavior of the dynamics, due to the nonlinear dependence, A = A(V{t)). Thus, 
the intricate aspect of the self-organized dynamics is tracing its algebraic connectivity 

over time through the self-propelled mechanism in which the nonlinear dynamics and algebraic 
connectivity are tied together. This issue is explored in the next section, dealing with local 
interactions and clustering. 



3. Local interactions and clustering 



In this section we consider the self-organized dynamics (1.1) of a "crowd" of N agents 



V — {pij^Li which does not interact globally: entries in their adjacency matrix may vanish, 
dij > 0. The dynamics is dictated by local interactions and its large time behavior leads to 
the formation of one or more clusters. 



^We ignore possible cases in which the self-organized dynamics may regain connectivity under "cluster dy- 
mics" , namely, agents separated into disconnected clusters and merging into each other at a later stage. 
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3.1. The formation of clusters. A cluster C is a connected subset of agents, {pi}ieC, which 
is separated from all other agents outside C, namely 

#1. dij for all i, j E C; and #2. — whenever i E C and j ^ C. 

The important feature of such clusters is their self-contained dynamics in the sense that 

j t Vi = J2^ij(Pj - Pi), ^2 a ij = 1 3 ^ E C. 

jeC jeC 

The dynamics of such self-contained clusters is covered by the concentration statements of 
global dynamics in section [2j In particular, if cluster C{t) remains connected and isolated for 
sufficiently long time, then its agents will tend to concentrate around a local consensus, 

Pi(*) ^Pc 5 , fora11 

The intricate aspect, however, is the last if statement: the evolution of agents in a cluster C 
may become influenced by non-C agents, and in particular, different clusters may merge over 
time. 

In the following, we fix our attention on the particular models for opinion and flocking dynam- 



ics, expressed in the unified framework (1.4), 



d N 1 

(3.1a) — pi = ^2,aij(pj -p^, dij = a^-(x) = — 0(|xf -x^-|); 

j=i 

Recall that pH^xin opinion dynamics, p^xin flocking dynamics, and a\ is the degree, 

{Ui — N : symmetric model, 

Gi — X^j/i 0(l x i ~~ non symmetric model. 

We assume that the influence function <j) is compactly supported 
(3.2) Supp{<f>(-)} = [0,R}. 

A cluster C = C(t) C {1, 2, . . . , N} is dictated by the finite diameter of the influence function 
(j) such that the following two properties hold: 

#1. max |xi(t) — x 7 -(t)| < R; and #2. min |x^(t) — x 7 -(t)| > R. 

i,jeC(t) ieC(t),j(£C(t) 

When the dynamics is global, R » [x(0)], then the whole crowd of agents can be considered 
as one connected cluster. Here we consider the local dynamics when R is small enough relative 



to the active diameter of the global dynamics: R < [x(0)] in the opinion dynamics (1.2), or 



R < Dqq in the flocking dynamics (1.3). The statements of global concentration towards a 
consensus state asserted in propositions 2.4| and[23]do not apply. Instead, the local dynamics 



of agents leads them to concentrate in one or several clusters — consult for example, figures 



1.1| and |3.3[ |5.2| below. Our primary interest is in the large time behavior of such clusters. 
The generic scenario is a crowd of agents which is partitioned into a collection of clusters, 
C/e, k = 1, . . . K, such that 

either \pi(t) - Pj(t) \ 0, if z, j |x^(t) - Xj(i)| < R 

or \xi(t)-Xj(t)\> R, if ieCkJeCi, k^L 
In this context, we raise the following two fundamental questions. 
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Question #1. Identify the class of initial configurations, V(0), which evolve into finitely 
many clusters, Ck, k = 1, . . . K. In particular, characterize the number of such clusters K for 
t » 1. 

Question #2. Assume that the initial configuration V(0)) is connected. Characterize the 
initial configuration V(0) which evolve into one cluster, K{t) = 1 for t ^> 1, namely, the 
question of emerging of consensus in the local dynamics. 

A complete answer to these questions should provide an extremely interesting insight into 
local processes of self-organized dynamics, with many applications. In the next two sections 
we provide partial answers to these questions. We begin with the first result which shows that 
if the solution of (3.1) has bounded time- variation then it must be partitioned into a collection 
of clusters. 

Proposition 3.1 (Formation of clusters). Let V(t) = {pk(t)}k be the solution of the opinion 
or flocking models ( |3.1| ) with compactly supported influence function Supp{(p(-)} = [0,i?) ; and 
assume it has a bounded time-variation 



/oo 
\Pi(s)\d 



s < oo. 

Then V{t) approaches a stationary state, p°° , which is partitioned into K clusters, {Ck\k=v 
such that {1, 2, . . . , N] = U% =1 C k and 

either Pi(t) pg°, for all i E Ck, 

(3.4) I 

k or |x°° - x°°| >R, if i g C k ,j €C £ , k + L 

Remark. Observe that if the solution decays fast enough — in particular, ifp(t) decays expo- 
nentially fast, \pi(t)\ < e -C(*-*o) j t > to > o ( as in the unconditional consensus and flocking 
of global interactions discussed in section [2]) , then it has a bounded time-variation. 

Proof Assumption ( |3.3[ ) implies 

\Pi(t 2 ) - Pi(h)\ < / \M s )\ ds < 1 for *2 > h > 1, 
Jtl 

hence each agent approaches its own stationary state, Pi(t) ? p°°. Therefore, since the 



expression on the right of ( |3.l[ ), 

p<(*) = Z>(M*) -Mt)\)(Pi(t) -Pj(t)), <n(t) = J^0(|xi(t) -x^t)!), 



3 



has a limit, it follows that lim^oo Pi (X) exists and by (3.3) it must be zero, Pi(t) —> 0. 
Take the scalar product of ( 3.1a| ) against p^ and sum, 

(3.5) ^2°i(Pi,Pi) = a^(j)ij(pj -p^ Pi) = ~^^2^ij\Pj ~Pi\ 2 - 

i ij ij 

The expression on the left tends to zero: p^ E Q(0) and G{ < N are uniformly bounded and 
Pi(t) — > 0; hence the limiting expression on the right yields 

(3.6) <K|x°° - xf |)|p- - p°°| 2 = 0, for all i,j < N. 

Thus, if |x?° — x.j°\ > i?, then agents i and j are in separate clusters. Otherwise, when they 
are in the same cluster, say z, j E Ck so that |x?° — x?°| < i?, then 0(|x?° — x?° |) > and by 
(3.6) they must share the same stationary state, p°° = p°° =: pg°, that is, (3.4) holds. □ 



HETEROPHILIOUS DYNAMICS ENHANCES CONSENSUS 



15 



We now turn our attention to the number of clusters, K. 



3.2. How many clusters? Note that if p°° = (pf , . . . , P^) T be a stationary state of (3.1) 
then p°° is an eigenvector associated with the nonlinear eigenvalue problem, 



A(x 00 )p 00 = p°°, 

corresponding to the eigenvalue \n(A(x.°°)) = 1. Actually, the number of stationary clusters 
can be directly computed from the multiplicity of leading spectral eigenvalues of \n(A(x.°°)). 

Proposition 3.2. Assume that the crowd of N agents {pi(t)}fL 1 is partitioned into K clusters, 
{1, 2, . . . , N} = U^^Cfc. Then, the number of clusters, K — K(t), equals the geometric 
multiplicity of AjvCA(x(t)) = 1, 

(3.7) K(t) = {#Ajv(A(x(t)) | A*(A(x(t)) = 1} . 

Pr oof. We include the rather standard argument for completeness. Suppose that the dynamics 



of (3.1) at time t consists of K — K(t) clusters, U^^C&. Define the vector r k = (rj, . . . , r^ T 



such that: 

k _ \ 1 if 3 € Ck 



We obtain 



y . 

3 otherwise 



dij. 

3 3^C k 

Using the fact that A is a stochastic matrix and that — if x^ and Xj are not in the same 
cluster, we deduce 

1 if z E C fc 
otherwise 



Y*k 

L i 5 



5^ a U = 

and therefore Ar k = r k . Thus, associated with each cluster there is an eigenvector r fc 
corresponding to Xn(A) = 1. To conclude the proof, we have to show that there are no other 
vectors r satisfying Ay — r. Indeed, assume that Ay — r, 

^CLijrj = n for any z., 

3 

Fix a cluster Ck- Then for any p E we have 

^ flp^rj = r p for any p e C k - 

Denote by r q the maximal entry of r^-'s on the left, corresponding to some q £ Ck'- since 
~ ^ with a p j > 0, we deduce that for any p £ Ck we have r p = a pj r j < Yl a pj r q — 
r q . Thus, the entries of r are constant on the cluster so that r oc Y k . □ 

3.3. Numerical simulations with local dynamics. We illustrate the emergence of clusters 
with one- and two-dimensional simulations of the opinions dynamics model (1.2b), 

(3.8) i*£_^ (xi _ Xi) , Xi (t)eK d 

ai ■ Z^k Pik 

The influence function, 0, was taken as the characteristic function of the interval [0, 1]: <j)(r) — 
X[o,i], and we use the Runge-Kutta method of order 4 with a time step of At = .05, for the 
time discretization of the system of ODEs (3.8). 
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Evolution in time of the opinions Xj 



* 6 




10 15 20 25 30 35 40 

time (t) 



Distribution of opinions Xj at t = 40 



Figure 3.1. The opinion model (|1.2b|) with M = 100 agents and <j) = X[o,i] (Left 
figure) and the histogram of the distribution of at t = 40 unit time (Right figure). 
We observe the formation of 4 clusters separated by a distance greater than 1. 



As a first example, we run a simulation of the one-dimensional opinion model, d = 1, subject 
to initial configuration of N = 100 agents uniformly distributed on the interval [0, 10]. In the 
figure [3TT] (Left) , we plot the evolution of the opinions x^(t) in time. We observe the formation 
of 4 clusters after 15 unit time. The histogram of the distribution of agents at the final time 
t — 40 (figure |3.1| right) shows that the distance between the clusters is greater than 1 as 
predicted by proposition |3.1| We also observe that the number of opinions contained in each 
cluster differs (respectively 35, 14, 31 and 20 agents). Indeed, the larger cluster at x « 2 with 
35 opinions is a merge between 3 branches (figure |3T| ) with one branch in the middle connecting 
the two external branches. When the two external branches finally connect at t « 8.5 (their 
distance is less than 1), we observe an abrupt change in the dynamics following by a merge of 
the 3 branches into a single cluster. 

To analyze the cluster formation, we also look at the evolution of the eigenvalues of the 
matrix of interaction A(x(t)) in (3.8), — 4>ij / J^k^ik- I n the- figure 3.2, we represent the 
evolution of the 8 first eigenvalues of the matrix A. From t = to t « 3.5, we observe that 
the 4 first eigenvalues converge to 1 which counts for the fact that only 4 clusters remains at 
this time. Then the matrix A(x(t)) remains constant in time from t w 3.5 to t « 8.5. At 
t « 8.5, two branches (see fig ure |3.1[ ) re-connect, and the two eigenvalues A5 and A6 equal 
zero. This confirms proposition 3.2 where the additional multiplicity of the spectral eigenvalue 
\jsr(A(-x(t)) = 1, indicates the formation of a new cluster. 

Next we turn to illustrate the dynamics of the two-dimensional, d = 2, opinion model 
( |1.2b ). With this aim, we run the model starting with an initial condition of N = 1000 agents 
distributed uniformly on the square [0, 10] x [0, 10]. We present, in figure [3T3| several snapshots 
of the simulations at different time (t = 0, 2, 4, 6, 12 and 30 unit time). As in the ID case, 
we first observe a fast transition to a cluster formation (from t = to t = 6). However at 
time t = 12, the dynamics does not have yet converged to a stationary state, we observe at 
the upper-left that three branches are at distance less than 1. This scenario is similar to the 
one observed in figure [3TT] with the apparition of 3 branches. At t = 30, the three clusters at 
the upper left have finally merged and the system has reached a stationary state: each cluster 
it at distance greater than 1 from each other. 
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- |A 2 

- |A 3 
- |A 4 

- |A 5 
|A 6 

- |A 7 

- |A 8 



10 15 

time (t) 



FIGURE 3.2. Absolute values of the eigenvalues of the matrix A (3.8) during the 



simulation given in figure |3.1| The number of eigenvalues equal to 1 corresponds to the 
number of clusters. 



4. K = 1: UNIFORM CONNECTIVITY IMPLIES CONSENSUS 



The emergence of a consensus in the opinion or flocking models (1.4) implies that the 
underlying graph associated with the dynamics must remain connected, namely |x^(t)— x^(t)| <C 
R at least for t ^> 1. Here we show that the converse is true: uniform connectivity implies 
consensus. Recall that the graph associated with (1.1), Qa '= (V,A(V)), is connected if every 
two agents Pi(t) and Pj(t) are connected through a path := {k\ = i < k<i < . . • < k r = j} 
of length rij < N. We measure the uniform connectivity by its "weakest link". 



Definition (Uniform connectivity). The self- organized dynamics (1.1) is connected if there 
exists ji{t) > such that for all paths Yij, 

(4.1) min a klM+1 {V{t)) > ^(t) > 0, for all ij. 

In particular, if /i(t) > \i > then we stay that V{t) is uniformly connected. 



Alternatively, uniform connectivity of (1.1 ) requires the existence of /i = ha > independent 
of time, such that 

{A N (P(t))) ij >n N >0. 

4.1. Consensus in local dynamics — symmetric models. We consider the symmetric 
dynamics ( |1.1[ ) with associated graph Qa •= (*PiA(V)). Fix the positions of any two agents 
Pi(t) and Pj(t) and their (shortest) connecting path Y^ or length r^. If we let diam(GA) '•— 
max^j r^ denote the diameter of the graph Qa, then 

I Pi - Pjl 2 < diam{G A ) ^ Ip**+i ~ P^J 2 > diam(g A ) < N. 
By uniform connectivity, fi < a^ +lj ^ along each path and hence 
( 4 - 2 ) &m(gJ PrP ^ ~ 5^ a k i+ ite\Pk l+1 -Pk t \ 2 <^a>ij\Pi-Pj\ 2 , 
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2 4 6 

Opinions at t 



10 



2 4 6 E 

Opinions at t = 6 



10 



2468 10 02468 

Opinions x^ at t = 12 Opinions x^ at t = 30 



FIGURE 3.3. Simulation of the opinion model ( |1.2b[ ) in 2D with M = 1000 agents 
and (j) = X[o,i]- The dynamics converges to a cluster formation (17 clusters) with each 
cluster separated by a distance greater than 1. 
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and summation over all pairs yields 



/A ^ l Pi ~ Pj 1 2 < N 2 y2aij\pi - pj\ 2 . 



diam(QA) 

ij ij 



Now we recall our notation := — (p): invoking (2.19) we find, 



|2 



wr x • (L A cj,q) . (1/2) E^^ji^-Pjl . /i 
(4.3) \2{La) mm — — - mm — — — — -==^ — . > 



Eq fc =o (q, q) p (1/2N) ■ | Pi - Pj f " Ndiam(g A ) ' 

Thus, the scaled connectivity factor \i j '(N diam(Q a)) < l^/N 2 serves as a lower bound for the 
Fiedler number associated with the symmetric dynamics of ( |1.1| ), |51j . Using theorem 2.6 
conclude the following. 



we 



Theorem 4.1 (Connectivity implies consensus: the symmetric case). Let V(t) = {pk(t)}k be 
the solution of a symmetric self- organized dynamics 

^Pi(t) = <* J^a«(^(*))(Pj(*) -P»(*))j a u = a ji- 

IfV(t) remains connected in time with "sufficiently strong" connectivity /j>a(v(s)) > 0; ^ 
approaches the consensus (p)(0) ; namely, 

V p( t) < exp ^ V p( o), Vj (t) := ^ ^ | Pi (t) - (p)(0)| 2 . 



In particular, ifV(t) remains uniformly connected in time, ( |4.1[ ) ; t/ien approaches an emerg- 
ing consensus, Pi(t) — >► p°° = (p)(0) a convergence rate, 

-a—t 

(4.4) V p(t) <e °^ 2 V p(0 ). 

It is important to notice that theorem |4.1| requires the intensity of connectivity to be suffi- 
ciently strong: connectivity alone, with a rapidly decaying /i(t), is not a sufficient for consensus 
as illustrated by the following. 

Counterexample. Consider the symmetric dynamics ( |1.2a[ ) with 5 agents, xi, . . . , X5, subject 
to initial configuration 

(4.5) xi (0) = -x 5 (0), x 2 (0) = -x 4 (0), x 3 (0) = 0, 

with (x4(0),X5(0)) to be specified below inside the box V := < X4 < 1 < X5 < |}. We 
fix the influence function <f>(r) = (1 + r) 2 (l — r) 2 X[o,i] 5 compactly supported on [0,1]; note 
that (f) f {fS) = 0'(1) = 0. By symmetry, the initial ordering in ( |4.5| ) is preserved in time. In 
particular, xs(t) = 0, and x 4 (t),X5(t) preserve the original ordering, \ < x 4 (£) < 1 < xs(t), 
the symmetric opinion dynamics Ql.2a ) (with a = 5 for simplicity), is reduced to 

/ 46 x x 4 = -0(|x 4 |)x 4 + 0(|x 5 -x 4 |)(x 5 -x 4 ) 

x 5 = 0(|x 4 - x 5 |)(x 4 - x 5 ) 

An equilibrium for the system is given by x 4 = X5 = 1. The eigenvalues of the linearized 
system at (1, 1) are Ai = and A2 = —2, therefore the equilibrium is unstable. By the 
unstable manifold theorem, there exists a curve, A4, tangent to the eigenvector V2 = (1, — 1) T 
at (1,1) such that any solution (x 4 (t), xs(t)) starting on M converges toward the unstable 
equilibrium (1,1) (see figure [XT] ) . Taking as initial condition (x 4 (0), xs(0)) onMflP yields 
a solution which converges toward 1. Since x 3 (£) = and < x 4 < 1, the configuration 
V{t) = {xi(t), . . . , X5(t)} stays connected in time (e.g. ji{t) := min^ <fii,i+i(t) > 0) but does not 
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converge to a consensus. The theorem |4.2| does not apply here because of the rapid decay in 
the strength of connectivity — /i(t) decays exponentially fast so that f°° /^(s)ds < oc, which 
is not sufficient to enforce consensus. 




Figure 4.1. Left: a solution of the symmetric model ( 1.2 ) that stays connected 
but does not converge to a consensus. Right: the initial configuration is picked 
from the unstable manifold A4 of the system (4.6). 



4.2. Consensus in local non-symmetric opinion dynamics. Next we turn to consider 
the question of consensus for the non-symmetric opinion model (1.2b). 

Theorem 4.2 (Connectivity implies consensus: non-symmetric opinion dynamics). LetV{t) = 
{Pk(t)}k be the solution of the non-symmetric opinion dynamics (1.2b) with compactly sup- 
ported influence function, Supp{(f)(-)} = [0, R), 

j k 

If V{t) remains uniformly connected in time in the sense that each pair of agents (z,j) is 
connected through a path Yij such tha^ 

min </>(\xk t ~ x ^ +1 |) > M > 0, for all ij, 

kieTij 

then it has bounded time-variation and consequently, V{t) approaches an emerging consensus, 
~Xi(t) — » x°° with a convergence rate, 

(4.7) |xi(t) - x°°| < e- am ^°)[x(0)], m = min^(r) > 0. 

Proof. We introduce the energy functional, 

£(t) :=a V$(|x,(t) -Xi(t)l), *(r) := f s<t>(s)ds, 
id Js=0 



2 Observe that here we measure connectivity in terms of the influence function <pij rather than the adjacency 
matrix, cnj as (4.1); the two are equivalent up to obvious scaling of the degree g%. 
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which is decreasing in time 
d 
dt 



(4.8) (t) = ( x j _ *i ? x j _ x *) = X] ^ ' x -?' ~ x * 



M h3 

= -2^(xi (xj - Xi)) = -2^a i |x i | 2 < 0. 

To upperbound the expression on the right of (4.8), sum ( |1.2b| ) against x^ to find, 

( 4 - 9 ) ^ X^'I X * ~ X ^'| 2 = -a^0ij(Xi - X^,Xi) = ^t7i(Xi,Xi) 



< ^y^ai\yLi\ 2 ^y^aj\i<Li\ 2 < N max | x^ ( 0) y ^ | Xj | 2 . 
We end up with the energy decay 

\ 2 

^ « \ 1 



(4.10) (t) < -\c?Cl ^^|x, - Xj fj 



iVmax^ |xi(0)| 



Hence, since 



J°° |X>^)i x *(*) -xj(t)\ 2 j dt < -J^m < oo, 

the sum 0r/(^)l x i(^) — x j(^)| 2 must become arbitrarily small at some point of time, namely, 
there exists to > such that 

(4-11) X>i(*0)l*i(*o) -x^o)| 2 < j^R 2 , 

and by uniform connectivity, consult ( |4.2[ ), 

(4.12) ^lxi(to) -x^to)! 2 < <f>k e ,k e+ M\xk t (to)-Xk t+1 (to)\ 2 < ^R 2 - 

Thus, the dynamics at time to concentrate so that its diameter, [x(to)] = niax^j |xj(to) — 
Xj(to)| < R/2, and since [x(-)] is non-increasing in time, [x(t)] < R/2 thereafter. Arguing 
along the lines of proposition |2.4| we conclude that there is an exponential time decay, 

Nctij > c/)(\xi(t) — Xj(t)|) > min <p(r) > min <p(r) = rn. t > to, 

K[x(t)] r<R/2 

and consensus follows from corollary (|2.3|). □ 



The decreasing energy functional £(t) can be used to estimate the first "arrival" time of 
concentration to. To this end, observe that: 

r\xj-*i\ r\xj-*i\ |x--x-| 2 
$(|x 7 - - xJ) = / s(/)(s)ds < M sds = M— M := max0(r). 

J 8=0 Js=0 2 

Using the assumption of uniform connectivity, there exists \i > and a path such that: 



|Xj — Xi| < — - <fiki,k i+1 I x ^ + i — x /cJ < ~~ / v ^,jl X j ~~ x ^ 



2 
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Combining the last two inequalities, we can upperbound the energy £: 



x MN 3 , 



Hence, (4.10) implies the Riccati equation 

d 



1 



dt m < 



which shows the energy decay 



£(t) 



< 



( 4/x 

\M7V 3 

1 



TV 6 



1 + 



Thus, the arrival time of concentration to (4.11) is at most of the order of 0(N 7 / fi 3 ). This 
bound on the first arrival time can be improved. 

We close this section by noting the lack of a consensus proof for our non-symmetric model 



of flocking dynamics (1.3b) is due to the lack of a proper decreasing energy functional. 



5. Heterophilous dynamics enhances consensus 

As we noted earlier, the large-time behavior of local models for self-organized dynamics 
depend on the details of the interactions, {a^-}, and in the particular case of local models 



(3.1), on the profile of the compactly supported influence function (j). Here we explore how the 



profile of (j) dictates cluster formation in the opinion dynamics model (1.2b). The numerical 



simulations presented in this section leads to the main conclusion that an increasing profile of 
(j) reduces the number of clusters {Ck\k=v ^ n particular, if the profile of (j) is increasing fast 
enough, then K = 1; thus, heterophilous dynamics enhances the emergence of consensus. 

In the following, we employ a compactly supported influence function <j) which is a simple 
step function, 



(5.1) 



r = < 



a 
b 




for r < j= 
for 4^ < r < 1 
for r > 1. 



The essential quantity here is the ration b/a which measures the balance between the influence 
of "far" and "close" neighbors (see figure |5T| ) . We initiate the opinion dynamics (1.2b) with 
random initial configuration {x^(0)}^. 



5.1. ID simulations. We begin with four simulations of the ID opinion dynamics (1.2b) 



subject to 100 opinions distributed uniformly on [0,10], the same initial configuration as in 
figure 3.1 To explore the impact of the influence step function (5.1) on the dynamics 



we 



used four different ratios of b/a = .1, 1, 2 and 10. As b/a increases, we reduce the influence of 
the closer neighbors and increase the influence of neighbors further away; thus, increasing b/a 



reflects the tendency to "bond with the other" . As observed in figure |5.2[ the increase in the 
ratio b/a — .1, 1, 2 and 10, reduces the corresponding number of limit clusters to K — 6, 4, 2, 
and for b/a — 10, the dynamics converged to a consensus, K = 1. The simulations of figure 



5.2 indicate that reducing the influence of closer neighbors and hence increasing the weight for 
the influence of neighbors further away, will favor increased connectivity and the emergence of 
consensus. 
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FIGURE 5.1. Influence functions <j) used in the simulations. The larger b/a is, the 
more heterophilous is the dynamics. 



To make a systematic analysis of the cluster formation dependence on the ratio b/a, we 
made several simulations with random initial conditions for a given ratio b/a. Then we make 
an average of the number of clusters, denoted by (S), at the end of each simulation (t — 100). 
To compute the number of clusters, we estimate the number of connected components of the 
matrix A ( |3.8[ ) using a depth- first search algorithm. As observed in figure |5.3| the number of 
clusters (S) decreases as b/a increases. Moreover, (S) approaches 1 when b/a approaches 10, 
implying that a consensus is likely to occur when b/a is large enough. 



5.2. Clusters and branches. We revisit the opinion model Ql.2b ) with an influence step 



function (5.1). As noted before, the increasing value of b/a increases the probability to reach a 



consensus. The simulations in figure 5.2 with b/a — 2 and with b/a — 10, show the apparition 
of branches, where subgroups of agents have converged to the same opinion yet, in contrast to 
clustering, these branches of opinions are still interacting with outsiders, which are in distance 
which is strictly less than R = 1. In particular, when b/a — 10, the distribution of opinions 
{x.i(t)}i, aggregate to form distinct branches seen in figure [5\2| at t ~ 5, one can identify in 
figure |5~H the formation of 10 branches which are separated by a distance of approximately .7 
spatial units. Since the distance between two such branches is always less than the diameter 
R = 1 of 0, these branches are not qualified as isolated clusters, as they continue to be 
influenced by "outsiders" from the nearby branches. Over time, these branches merge into 
each other before they emerge into one final cluster, the consensus, at t ~ 33. Thus, the 
decisive factor in the consensus dynamics is not the number of branches but their large time 



connected components. Indeed, figure 5.2 with b/a — 10, shows that the agents in the different 
branches remain in the same connected component at distance ~ .7, corresponding to the 
discontinuity of </>(•), which experience a jump from .1 to 1 at l/y/2 « .7. 

To illustrate the apparition of the distance l/y/2 between two nearest branches, we repeat 
the simulations, this time with special initial configurations where all the opinions are uniformly 



spaced with |x^+i — x^| = with < < 1. As we observe in figure 5.5, the agents {x^ 
readjust their "opinion" such that the distance between nearest neighbors |x^ + i— x^| approaches 
l/y/2 as t > 1. 

5.3. 2D simulations. We made several 2D simulations with different influence functions (j). 
As a first illustration, we made a 2D simulation with the same initial configuration used in 



figure 3.3, but this time we used the influence step function (j) in (5.1) with b/a — 10. In figure 



5.6, one can observe a concentration phenomenon (from t = to t = 2.5) — the opinions 
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55^ 



with b/a = .1 



t 4> 



5 10 15 20 25 30 35 40 

time (t) 



X 6 



o 




(j) with b/a = 1 



t 



o 1 



5 10 15 20 25 30 35 40 

time (t) 



cj) with b/a = 2 




10 15 20 25 30 35 40 

time (t) 



with b/a = 10 

1 J 




FIGURE 5.2. Simulation of the opinion dynamics model with different interacting 
function (j). When the influence of close neighbors is reduced (i.e. b/a large), the 
number of cluster decreases. For b/a = 10, the dynamics converges to a consensus. 



aggregate into 5 final clusters, compared with the 17 clusters observed in figure [373] with the 
influence function <j) — X[o,i]- Thus, as in the ID case, a more heterophilous influence function 
increases the clustering effect. 

We also estimate the average number of clusters (S) depending on the ratio b/a. As observed 
in figure [pVT} (S) is a decreasing function of b/a and once again the decay of (S) as a function 
of b/a E [0, 10] is logarithmic. 



6. Heterophilous dynamics - nearest neighbor model 

Careful observations of startling flocks led the Rome group jT5l [T6l [T7] to the fundamental 
conclusion that their dynamics is driven by local interaction with fixed number of nearest 
neighbors. This motivates our study of nearest neighbor models for opinion dynamics which 
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ratio 6/a 
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FIGURE 5.3. Average number of clusters (S) depending on the ratio b/a (Left figure). 
The larger b/a is, the fewer the number of clusters. The decay is logarithmic on [1, 10] 
(Right figure). For each value of b/a, we run 100 simulations to estimate the mean 
number of clusters (S). Simulations are run with At = .05 and a final time equals to 
t = 100 unit time. 

Distribution of opinions x.; at t = 5 



FIGURE 5.4. The distribution of {x^ in the simulation of figure [5^2] with b/a — 10 
at time t = 5. The distance between two picks of density is around l/y/2 ~ .7 space 
unit. This distance corresponds to the discontinuity of the function 4>(r). 



take the form 
(6.1a) 



d v-^ 

dt ^ a, 

h- \j-A<q} 



^( Xi 

where the degree Oi is given by one of two forms, depending on the symmetric and non- 



symmetric version of the opinion dynamics in (1.2) 

the symmetric case : 



(6.1b) 



(Ji 



the nonsymmetric case : G{ 



1 

2q 

{j-\j-i\<q} 



Thus, each agent i is assumed to interact only with its 2q agents i — q, . . . ,i + q. Typically, 
q is small (the observation in [T5| [T7] report on six to seven active nearest neighbors). 



26 



SEBASTIEN MOTSCH AND EITAN TADMOR 




FIGURE 5.5. Initial condition with an equi-repartition of {x^: |x^ + i — x^| = .9 (Left 
figure) and |x i+ i — x^| = .6 (Right figure). Nearest neighbor readjust their distance to 
l/\/2 « .7 unit space, we observe a concentration of the trajectories in the left figure 
and a spread of the trajectories in the right figure. 



We analyze the connectivity of the particular case of two nearest neighbors, q = 1. Here we 
prove that such local models preserve connectivity and hence converge to a consensus provided 
the influence function <j) is increasing. This result supports our findings in section [5] that 
heterophilous dynamics is an efficient strategy to reach a consensus. 

6.1. Nearest neighbor with global influence function. We begin by noting that the 
different approaches for consensus of global models apply in the present framework of local 
nearest neighbor models (6.1). For example, consider the non-symmetric nearest neighbor 
model 

(6.2) Gi< Jt* i = a 5^ ^j( x j ~ x *)' a i = ^2 far 

{j--\j-i\<q} {j--\j-i\<q} 
It admits an energy functional, 

£(*):= a $(| Xj (i)-Xi(i)|), $( r ) := s<f>(s)ds, 

{i j\<q} ' S ~° 

which is decreasing in time, £(t) < £(0) and we conclude 



Theorem 6.1. (Global connectivity) Consider the nearest neighbor model (6.2) with an influ- 
ence function 0, Supp{(j)(-)} = [0, R) and assume a<f>(R) > £(0)- Then min^_j^ <q 4>ij(t) > 
where := min r<j R0(r). Hence, the nearest neighbor dynamics (6.2) remains connected and 
consensus follows. 

Proof. Since £ is decreasing in time, 

a$(|xi(t) -Xj(t)|) < £(0) < a$(R) for any \i - j\ < q, 

and since <&(r) = J r sc/>(s)ds is an increasing function, |xj(t) — Xj(t)| < i?, hence faj > 
and consensus follows. □ 

We note, however, that since < <j> < 1, then $( r) h as a quadratic bounds, m^r 2 < 
2$(r) < r 2 , and hence the assumption made in theorem |6. l| implies 

a$(R) > £(0) ~> R 2 >m OQ ^ |x* - Xj | 2 . 

\i-j\<q 
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FIGURE 5.6. The heterophilous effect: diminishing the influence of close neigh- 
bors relative to those further away, increases the clustering effect. 2D simulation 



of the opinion model (1.2b) with M = 1000 agents using a step influence function 
<j) = .1 X[q i/y/2] + X[i/y2 i] leads to 5 clusters which remains at the end of the simula- 



tion. This should be compared with 17 clusters with cf) = X[o,i] ( see figure 3.3). 
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ratio b/a ratio b/a 



FIGURE 5.7. Average number of clusters (S) depending on the ratio b/a in 2D (Left 
figure). As in ID case, the larger b/a is, the fewer the number of clusters, and the decay 
is logarithmic on [1, 10] (Right figure). For each value of b/a, we made 100 simulations 
to estimate the mean number of clusters (S). Simulations were made with At = .05 
and were recorded at the final time t = 100. 



Namely, the support of (j) should be sufficiently large to cover a globally connected path in 
phase space. 

6.2. Two nearest neighbor dynamics. In this section we prove uniform connectivity and 

hence convergence to a consensus of a symmetric two nearest neighbor model, (6.1), 

(6.3) 

d _ at , . ,\ _ J 0, i = 0,N 

We assume that the initial configuration of agents can be enumerated such that {x^(0)}i is 
connected 

(6.4) max|x z+1 (0) -x,(0)| < R, Supp{(/)(-)} = [0,12). 

i 

The configuration of such "purely" local interactions applies to the one-dimensional setup 
where each agent is initially connected to its left and right neighbors; we emphasize that these 
configurations are not necessarily restricted to the one dimensional setup. 

Forward differencing of (6.3) implies that A i+ i := A i+ i (t) := x^ + i(t) — x^(t) satisfy 

^A- + i = |(^ + |(x, +2 - Xi+i) + (Xi - Xi+i) - (Xi+i - X,) - (Xi_i - x, 



= ^(^sA^s - 2k <+ iA <+ ! + ^_|A-_i), i = l,2,...,JV. 



2 



The missing A's for i = \ and z = N + i are defined as A i = A 7V+ i = 0. Let A^ + i denote 
the maximal difference, |A p+ i | = max^ |A i+ i | measured in the ^2-norm. Then 



2 " 1 2 

1 d . A l9 Oi 



|2 

■p+41 • 



< f (« p+ | -2« p+ i+ N _i)|A 
Now, if is non- decreasing influence function, then 

\\ + ±\>\A l+h \ - 2« p+ i=2^(|A p+| |)>0(|A p _ | |) + ^(|A p+ 
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and hence |A + i(£)| 
following theorem. 



maxj cf)(\xi + i(t) — Xi(t)|) < maxj </>(|xj_|_i(0) — Xj(0)|). We deduce the 



Theorem 6.2. Consider the nearest neighbor dynamics (6.3) subject to initial configuration, 
x(0) which is connected, (6.4), 

max|x i+ i(0) -Xi(0)| < R, Supp{(f>{-)} = [0,R). 



Assume that the influence function, (f>, is non- decreasing. Then the dynamics (6.3) remains 
connected and converges to a consensus, x°° = (x)(0), 

2(p(0)t\ 



(6.5) 



E 



Xi(*) 



(x)(0)| 2 <exp 



N ) 



| Xl (0)-(x)(0)| 2 , (x>:=^E 



X 7 ; 



It is important to notice that theorem |6.2| requires an non-decreasing influence function. 
Indeed, the steeper the increase of <j) is, the better the connectivity is. This is concrete ram- 
ification of our main statement that heterophilious dynamics enhances consensus. Note that 
a two-nearest neighbor dynamics driven by a decreasing <j) will not guarantee consensus as 
illustrated by the following. 

Counterexample. We revisit the counterexample in section 4J_, of five agents symmetrically 
distributed around xs(t) = with \ < x 4 (t) < 1 < xs(t) < |, governed by 

x 4 = -0(|x 4 |)x 4 + 0(|x 5 - x 4 |)(x 5 - x 4 ), 

X 5 = 0(|X4 - X 5 |)(X 4 - X 5 ), 

with a compactly supported influence function <j>{r) — (1 — r) 2 (l + 2 X[o,i]- Observe that 
this configuration amounts to a two-nearest neighbor dynamics. Its concentration into three 
separate clusters { — 1, 0, 1} shown in figure 4.1, requires a rapidly decreasing influence function 
(to be precise — (j)(r)r I for r ~ 1), which is not covered by the two-nearest neighbors' 



heterophilious dynamics sought in theorem 6.2 



Proof. The adjacency matrix associated with (6.3), x = a (Ax — x) is given by the tridiagonal 
matrix A = {a^ }, given by 



K i+j , K i+j 
Z 2 2 



0(|Xi - Xj 



\ % -3\ 



1, 



The corresponding Laplacian associated with A is given by 



(6.6) 



La ~ — 

A 2 



K3 
2 



ft3 
2 

-«2 



K5 
2 



—^3 

2 



K 7 
2 



2 



Since |x^ + i(t) — x^(t)| < i?, the off-diagonal entries = 0(|x^ + i — x^|) > and hence, the 
graph Qa — ({x(t)} 3 A(x(t))) remains connected with /i = min^, i and diam(GA) — N. By 
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(Ob we find 



ML a) > 



mini 0(0) 



AT2 



AT2 



Using theorem 2.6 (see ( 2.21[ >) we end up with 



a^(0)t 

£ I^W " < x )(°)| 2 S e"^V2- £ | X . (0) _ (x)(0) | 2) 



which concludes the proof. 



□ 



Remark. The worst case scenario for the decaying of the \xi{t) — (x)\ is to have many opinions 
Xi concentrate at two extreme values with just one path of opinion connecting the two extremes 
(see figure 6.1). 



b - 



V2 



1 



FIGURE 6.1. The worst case scenario for the decaying of the norm of the vector A: 
the formation is connected but there are two large groups with extreme values. 



7. Self-alignment dynamics with discrete time steps 

Models for opinion dynamics were originally introduced as a discrete algorithms. In this 
section we therefore extend our results on the semi-discrete continuous opinion dynamics (1.2b) 
to the fully discrete case, 

[ ' At a Ej <t>ij 

In particular, for a = 1/ At we find that x™ = x^(nAt) satisfies the Krause model [H El [43] 
(7.2) x^ 1 = </hj = m ~ x?|). 



In the following, we study the properties of the discrete dynamics (7.2). 



HETEROPHILIOUS DYNAMICS ENHANCES CONSENSUS 



31 



7.1. Consensus with global interactions. Many results of the continuous dynamics (1.2b) 
can be translated to the discrete dynamics (7.2). For example, the convex hull of the opinions 
(2.3) is still decreasing in time: 

fi(ra + l) C Q(n). 



The discrete dynamics (7.2) will also converge to a consensus if initially all agents interact 



with each other. More precisely, arguing along the lines of proposition |2.4| gives the following 
result. 

Theorem 7.1. Assume that m = min rG [o 5 [ x (o)]] <j>(r) > 0. Then, the diameter of the discrete 
dynamics (7.2) satisfies 

(7.3) [x n ] < (l-m) n [x°] ™ 0. 

and convergence to a consensus, x™ n —^ ) x °° £ O(0) follows. 



Proof. Using the contraction estimate (2.7) followed by the bound tp A > max<9 9 • \{9) yield 

[x n+1 ]<(i-VA)[x n ]<(i-o-Ke,t n ))[* n ], » - ^ 



Fix 9 = m/N then A(9) includes all agents, \(9,t n ) = TV, and we conclude 

[x n+1 ] < (1-m) [x n ], 



which proves (7.3). 



□ 



7.2. Clustering with local interactions. As in the continuous dynamics, we would like to 
investigate the behavior of the discrete dynamics (7.2) with local interactions; in particular, 
we are interested in the formation of clusters. Our aim is to reproduce the discrete analog of 
proposition |3.1| 



Proposition 7.2. Let V n — {^}k be the solution of the discrete opinion dynamics (7.2) 
with compactly supported influence function Supp{<j){')} = [0, R). Assume that it approaches a 
steady state fast enough so that 



(7.4) 



EE 



x n+l_ x n|«™ a 



Then {x n } approaches a stationary state, x°°, which is partitioned into clusters, {Ck}k> such 
that {1, 2, . . . , N} = U% =1 C k and 



(7.5) 



for all i E 



Proof. By assumption (7.4) 



ri2 — 1 



ixr-*ri< E i x 



n+ 1 v n 



xf| < 1, for n 2 > ni > 1, 



n=n\ 



and hence x n approach a limit, x™ - 
the following form: 

j>,(x? +i -x?)= x;^-(x7-x?). 

3 3 



x?°. The discrete dynamics (7.2) can be written in 
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Taking the scalar product against xf, summing in i and using the symmetry of faj yields 



2 

u u u 

Since <^-, x™ and by assumption, the tail J2m l x r +1 — x ?l are bounded, we conclude that the 
sum on the right converges to zero 

0y|x? - x?| 2 ^ 0(xf - xfDIxf - xf| 2 = 0. 

Hence, either x?° and x?° are in separate clusters, |x^° — x?°| > R or else, they are in the 
limiting point of the same cluster, say i,j E Q so that x?° = x?°. □ 

We now turn our attention to the convergence toward consensus for the discrete dynamics 
(7.2). As for the continuous dynamics (1.2), there exists a Lyapunov functional energy for 
the dynamics under the additional assumption that the influence function <j) is non-increasing. 
Consequently, we deduce the analog of theorem |4.2| for the discrete dynamics. 



Theorem 7.3. Let V n — {^}k be the solution of the discrete opinion dynamics (7.2) with 
non-increasing, compactly supported influence function Supp{<j){')} = [0,i?). If V n remains 
uniformly connected for any n, then V n converges to a consensus. 

Proof. First, we prove that the energy functional £ n is also a Lyapunov function for the discrete 
dynamics: 



(7.6) £":=2>(|x?-x?|), *(r) = J\<f>(s)ds. 



Introducing (f(r 2 ) = $(r), we have (p(r) — s<p(s)ds = \ J Q r <p(- s /y)dy. By assumption (f) is 
non-increasing, thus ip is concave-down. Therefore, 



£ n + l_ £ n = ^ (p (| x «+l_ x n + l|2 ) _ ¥ , ( | x n_ x n|2 ) 



2 



Using |a| 2 — |b| 2 = (a — b , a + b), we deduce: 

£ n+1 -£ n < ^<^A t x?-A t x?,x? +1 -x? +1 +x?-x?) 



2 



since <fiij = <j)ji. Writing x™ +1 = x ™ + ^x™, we obtain: 

< ^«^<A t x?,2(x?-x?) + A t x?-A t x?}. 

Combining with the equality: 

(7-7) XI^ A ' x ? = E^( x i- x ?)- 
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we conclude 



n\2 
i I 5 



where we use once again the symmetry of the coefficients <pij. Thus, £ n is decaying. 

Now, we would like to combine the decay of £ n and the strong connectivity of V n . Noting 



Gi — ^2 j fiiji the equality (7.7) yields: 



Thus, 



iVmax^ \x{ (0)| 



Summing in n, we deduce that the sum ~~ x ?l becomes arbitrarily small. To 



conclude, we proceed as in the proof of theorem 4.2 



□ 



7.3. Numerical simulations of discrete dynamics. In this section we illustrate the differ- 
ence between the continuous opinion model (1.2b) and its discrete version (7.2). To this end, 



we run in parallel numerical simulations of the discrete and continuous model subject to the 
same initial conditions. 

First, we run a simulation with an influence function (j) = X[o,i] (figure 7.1). Discrete and 
continuous dynamics are very similar, except that there are three branches m the continuous 
dynamics which are not present in the discrete dynamics. For this reason, at the end of the 
simulation, we count 4 clusters in the discrete dynamics and only 3 in the continuous version. 

C[i/V2,i] with b / a 



Next we use the influence function (5.1) 



^[0,1/^] 



10. Here, the 

discrete and continuous dynamics give very different results shown in figure |7.2[ As we have 



seen previously, the continuous dynamics converges to a distribution with uniformly spaced 
clusters and then reach a consensus. In contrast, the discrete dynamics does not stabilize. 
Order between the opinions {x^ is no longer preserved, trajectories do cross. Even though 
the total number of clusters has been diminished with b/a = 10 (from 4 to 3 clusters), the 
effect of the ratio b/a on the clustering formation is less pronounced in the discrete dynamics. 



8. Mean-field limits: self-organized hydrodynamics 

When the number of agents N is large, it is convenient to describe the evolution of the 
resulting large dynamical systems as mean-field equation. We limit ourselves to a few classic 
general references on this topic [58l EH1 [33] , and a few recent references in the context of opinion 
hydrodynamics [62j HQ], and in flocking hydrodynamics fM [T^ [T^ [49| [52| [42] . 

8.1. Opinion hydrodynamics. To derive the mean- field limit of the opinion dynamics model 

empirical distribution p(£, 

1 N 

^' X ) := ArZX-(t)(x), 



(1.2b), we introduce the so-called empirical distribution p(£, x): 

N 

N 



3=1 
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Discrete model, b/a = 1 Continuous model, b/a = 1 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 



time (t) time (i) 

FIGURE 7.1. Simulations of the discrete (Left figure) and continuous dynamics 
(Right figure) with </> = X[o,i] starting with the same initial condition. Although 
the two simulations are very similar, the discrete dynamics yields 4 clusters whereas 
the continuous dynamics gives 3. We use a time discretization of At = .05 to simulate 
the continuous dynamics. 



Discrete model, b/a = 10 Continuous model, b/a = 10 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 



time (t) time (i) 

FIGURE 7.2. Simulations of the discrete (Left figure) and continuous dynamics 
(Right figure) with </> = •^- X\o l i / y /2] + X[i/V2 1] starting with the same initial condi- 
tion. In contrast with figure |77l| the two models produce very different output. There 
is no uniformly spaced formation in the discrete model, we only observe cluster forma- 
tion. 



where 8 is a Dirac mass and {xj(t)}j is the solution of the consensus model (1.2b). Expressed 
in terms of this empirical distribution, the non-symmetric model (1.2b) (with a = 1) reads, 



J y 0(|y- x il)(y- x i)p(^y) rf y _ 0(|y|)y *p)(xO 



/ y <K|y- x ;l)p(^y) d y 



(<K|y|)*p)(x.) 
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This equation describes the characteristics of the density p. Indeed, integrating p against a 
test function f yields^] 



dt 



(p, <p) 




\ 1 N 

<p(xj(t)) = -^(x^Vx^x^))). 



Using the expression (8.1), we deduce: 



it (t>-v) = jv^\ 4>(\y\)>p(*j) 



0(|y|)y *p 
<K|y|) *p ' 



0(|y|) *p 

Thus, p = x) satisfies a continuum transport equation, 



3.2a) 



^ + v x -(H = o with u ( x)= ; (|y _ x|)p(y)dy • 



This is the hydrodynamic description of the agent-based opinion model ( 1.2b| ). Similarly, the 
opinion hydrodynamics of the corresponding symmetric model (1.2a) (with a = 1) amounts to 
the aggregation model [TO] 



3.2b) 



dtp + V x • (pu) = with u(x) = V$ * p, $'(r 



roir). 



We note that the transport equations (8.2) are non-linear due to the dependence of the velocity 
field u = u(p). Main features of the particle description for opinion dynamics (1.2) carry over 
the hydrodynamic model ( |8.2[ ). Thus, for example, the symmetric mode l (|8.2b ) preserve the 
center of mass, (J x xp(t, x) dx) = where the non-symmetric model (8.2a) does not. We 



distinguish between the two cases of global and local interactions. 

The existence of regular solutions of the symmetric aggregation model (8.2b) for bounded 
decreasing 0's such that \(j)'{r)r\ < <p(r) was proved in [3]. This holds independently whether <j) 
is global or not. Moreover, if the kernel <j) is globally supported, then one can argue along the 
lines of the underlying agent-based model (8.1), to prove convergence of the hydrodynamics 
toward a consensus, that is, p(t, x) converges to a single point asymptotically in time. If <j) is 
compactly supported, however, then the velocity field u need not be continuous with respect 
to p due to the singularity when J y 0(|y — x|) p(y) dy = 0. Then, existence and uniqueness of 
solution of the non-symmetric model ( |8.2a| ) cannot be obtained through a standard Picard's 
iteration argument. The large time behavior of the dynamics in this local setup is completely 
open. As in the agent-based dynamics, the generic solution p(t, x) is expected to concentrate 
in a finitely many clusters, or "islands"; in particular, under appropriate assumption on the 
persistence of connectivity among these islands, one may expect a consensus. Preliminary 
simulations show that cluster formation tends to persist for the hydrodynamic model, but 
analytical justification remains open. 



8.2. Flocking hydrodynamics. We study the second-order flocking models (1.3) in terms of 
the empirical distribution f N (t, x, v) := X^=i ^xj(t)( x ) ® ^vj(t)( v )? where c5 x ® S v is the usual 



(•, •) denotes the duality bracket between distributions and test functions 
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Dirac mass on th e phase space i 1 * x Rf Consider the non-symmetric particle model system 
for flocking (1.3b): expressed in terms of f N , it reads 

d Xi dvi N J y ,w <Kly- x l) (w-vj) /(y, w) dydw 

dt dt J y ^(|y-x|)/(y,w)dydw 

which leads to Liouville's equation, 

(8-3) a t / + v-V x / + V v -(F[/]/) = 0. 

Integrating the empirical distribution in the velocity variable v yields the hydrodynamic 
description of flocking, expressed in terms of the density and momentum distributions of par- 
ticles, 

p(*> x ) = J /(*> x > v ) dv (corresponding to ^ 5 Xj (t)(x)) , 

p(t, x)u(t, x) = y v/(t, x, v) dv (corresponding to ^ v i (t)5 Xj . (t) (x 



Integrating the kinetic equation ( |8.3[ ) against the first moments (l,v) yields the system, cf., 



(8.4a) d t p + V x • (pu) = 

(8.4b) <9 t (pu) + V x • (pu ® u + P) = ap(u - u). 

The expression on the right reflects alignment: the tendency of agents with velocity u to relax 
towards the local average velocity, u(x), dictated by the normalized influence function a(x,y), 



i.4c) 



u(x) := / a(x,y)p(y)u(y)dy, / a(x, y)p(y) dy = 1. 



This includes in particular, the hydrodynamic description of the symmetric and non-symmetric 
flocking models, given respectively by 

f <Kly 



«(x,y) 



x 



C-S model (1.3a) 



<K|y-x| 



f 0(|y-x|)p(y)dy 



non-symmtric model (1.3b) 



The system (8.4) is not closed since the equation for pu ( 8.4b| ) does depend on the third 
moment of / which is encoded in the pressure term P := J v (v — u) ® (v — u)/(t 5 x, v) dv. 
If we neglect the pressure (in other words, assume a monophase distribution, /(t,x,v) = 
p(£, x) 5 u (t 5 x)( v ) so that P = O)? then the flocking hydrodynamics (8.4) is reduced to the closed 
system 

" dtp + V x • (pu) = 0, 

1.5) ; 

d t u + (u • V x )u a(u - u). 



The question of an emerging flock in ( |8.5[ ) follows along the lines of our discussion on the 
underlying agent-based models (1.3). The case of a global influence function is rather well- 



understood: in particular, regularity of the one-dimensional "incompressible" case, p = 1, 
depends on initial critical threshold [56l 26] . Flocking hydrodynamics governed by locally 
supported influence function requires a more intricate analysis, due to the realistic presence of 
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vacuum, [59]. The hydrodynamic description of self-organized dynamics give rise to systems 
like (8.5) which involve nonlocal means. Questions of regularity and quantitative behavior of 
such systems provide a rich source for future studies. 



References 

[1] I. Aoki, A simulation study on the schooling mechanism in fish Bull. Japanese Society of Scientific Fisheries, 
48(8), 1081-1088, 1982. 

[2] M. Ballerini, N. Cabibbo, R. Candelier, A. Cavagna, E. Cisbani, I. Giardina, V. Lecomte, A. Orlandi, G. 

Parisi, A. Procaccini, M. Viale and V. Zdravkovic, Interaction ruling animal collective behavior depends 

on topological rather than metric distance PNAS, 105(4), 1232-1237, 2008. 
[3] A. Bertozzi, J. Carrillo and T. Laurent, Blow-up in multidimensional aggregation equations with mildly 

singular interaction kernels Nonlinearity 22 (2009) 683-710. 
[4] V. Blondel, J. M Hendricks, A. Olshevsky, and J. Tsitsiklis, Convergence in multiagent coordination, 

consensus, and nocking. In IEEE Conference on Decision and Control, volume 44, page 2996, 2005. 
[5] V. D. Blondel, J. M. Hendricks, and J. N. Tsitsiklis, On Krause's multi-agent consensus model with state- 
dependent connectivity. Automatic Control, IEEE Transactions on, 54(11) : 2586-2597, 2009. 
[6] S. Boyd, A. Ghosh, B. Prabhakar, and D.Shah, Randomized gossip algorithms. IEEE Trans. Inform. 

Theory, 52, 2508-2530, 2006 
[7] A. M. Bruckstein, N. Cohen, A. Efrat, Ants, crickets and frogs in cyclic Pursuit CIS report #9105, Center 

for Intelligent Systems, Technion Israel Inst, of Tech., 1991. 
[8] S. Camazine, J. L. Deneubourg, N. R Franks, J. Sneyd, G. Theraulaz, and E. Bonabeau, Self-organization 

in biological systems. Princeton University Press; Princeton, NJ: 2001, 2001. 
[9] J. A. Canizo, J. A. Carrillo and J. Rosado, Collective behavior of animals: Swarming and complex patterns. 

2009. 

[10] C. Canuto, F. Fagnani, and P. Tilli, An Eulerian approach to the analysis of Krause's consensus models 

SIAM J. Control Optim., 50(1), 243-265. 2012. 
[11] R. Carli, F. Fagnani, A. Speranzon, and S. Zampieri, Communication constraints in the average consensus 

problem. Automatica, 44, 671-684, 2008. 
[12] J. A. Carillo, M. D'Orsogna, V. Panferov, Double milling in self-propelled swarms from kinetic theory. 

Kinet Relat Models, 2:363-378, 2009. 
[13] J. A. Carrillo, M. Fornasier, J. Rosado, and G. Toscani, Asymptotic flocking dynamics for the kinetic 

Cucker-Smale model. SIAM J. Math. Anal, 42:218-236, 2010. 
[14] J. A. Carrillo, M. Fornasier, G. Toscani, and F. Vecil, Particle, Kinetic, and Hydrodynamic Models of 

Swarming, in Naldi, G., Pareschi, L., Toscani, G. (eds.) Mathematical Modeling of Collective Behavior in 

Socio- Economic and Life Sciences, Series: Modelling and Simulation in Science and Technology, Birkhauser, 

(2010), 297-336. 

[15] A. Cavagna, I. Giardina, A. Orlandi, G. Parisi, A. Procac-cini, M. Viale and V. Zdravkovic, The starflag 
handbook on collective animal behaviour. 1: Empirical methods. Animal Behaviour 76:217-236, 2008. 

[16] A. Cavagna, I. Giardina, A. Orlandi, G. Parisi and A. Procaccini, The starflag handbook on collective 
animal behaviour. 2: Three-dimensional analysis. Animal Behaviour 76:237-248, 2008. 

[17] A. Cavagna, A. Cimarelli, I. Giardina, G. Parisi, R. Santagati, F. Stefanini and M. Viale, Scale-free 
correlations in starling flocks. Proc. Nat. Academy Sci., U.S.A. 107:11865-11870, 2010. 

[18] C. Cercignani, R. Illner, M. Pulvirenti, The mathematical theory of dilute gases. Springer series in Applied 
Mathematical Sciences, 106, Springer- Verlag, 1994. 

[19] D. Comaniciu and P. Meer, Mean Shift: A Robust Approach Toward Feature Space Analysis. IEEE Trans- 
actions on Pattern Analysis and Machine Intelligence archive Volume 24(5) Pages 603 - 619, 2002. 

[20] J. Cortes, S. Martinez, and F. Bullo, Robust rendezvous for mobile autonomous agents via proximity graphs 
in arbitrary dimensions. IEEE Trans. Automat. Control, 51, 1289-1298, 2006. 

[21] I.D. Couzin, J. Krause, N.R. Franks, and S.A. Levin, Effective leadership and decision-making in animal 
groups on the move. Nature, 433(7025) :513-516, 2005. 

[22] I. D. Couzin, J. Krause, R. James, G. D. Ruxton and N. R. Franks, Collective Memory and Spatial Sorting 
in Animal Groups J. of Theoretical Biology, 218(1), 1-11, 2002. 

[23] F. Cucker and S. Smale. Emergent behavior in flocks. IEEE Transactions on automatic control, 52(5) :852, 
2007. 



38 



SEBASTIEN MOTSCH AND EITAN TADMOR 



[24] F. Cucker and S. Smale, On the mathematics of emergence. Japanese Journal of Mathematics, 2(l):197-227, 
2007. 

F. Cucker, S. Smale and D.X. Zhou, Modeling language evolution. Found. Comput. Math., 4, 315-343, 2006. 
P. Degond and S. Motsch, Continuum limit of self-driven particles with orientation interaction Math. Models 
Methods Appl. Sci., 18(1):1193-1215, 2008. 

D'Orsogna, M. R., Chuang, Y.-L., Bertozzi, A. L. and Chayes, L., Self-propelled particles with soft-core 
interactions, patterns, stability, and collapse. Phys. Rev. Lett. 96, 104-302 (2006). 

B. Duering, P. Markowich, J.F. Pietschmann, and M.T. Wolfram, Boltzmann and Fokker-Planck equa- 
tions modelling opinion formation in the presence of strong leaders. Proceedings of the Royal Society A: 
Mathematical, Physical and Engineering Science, 465(21 12) :3687, 2009. 

L. Edelstein-Keshet, Mathematical models of swarming and social aggregation, International Symposium 
on Nonlinear Theory and its Applications, (NOLTA 2001) Miyagi, Japan, 2001. 

R. Eftimie, Hyperbolic and kinetic models for self-organized biological aggregations and movement: a brief 
review. J Math Biol. 65(l):35-75, 2012. 

F. Fagnani and S. Zampieri, Randomized consensus algorithms over large scale networks. IEEE J. Selected 
Areas of Communications , 26, 634-649, 2008. 

M. Fiedler, Algebraic connectivity of graphs, Czech. Math. J. 23(98), 1973, pp. 298-305. 
F. Golse, The mean- field limit for the dynamics of large particle systems. Journees Equations aux derives 
partielles, 9:1-47, 2003. 

V. Grimm and S. F Railsback, Individual-based modeling and ecology. Princeton Univ Pr, 2005. 
S. Y Ha and J. G Liu, A simple proof of the Cucker- Smale flocking dynamics and mean- field limit. Com- 
munications in Mathematical Sciences, 7(2):297-325, 2009. 

S. Y. Ha and E. Tadmor, From particle to kinetic and hydrodynamic descriptions of flocking. Kinetic and 
Related Models, l(3):415-435, 2008. 

R. Hegselmann and U. Krause, Opinion dynamics and bounded confidence: models, analysis and simulation. 
Journal of Artificial Societies and Social Simulation, 5(3), 2002. 

D. Helbing, Traffic and related self-driven many particle systems. Reviews of Modern Physics, 73, 1067-1141, 
2001. 

C. K Hemelrijk and H. Hildenbrandt, Self-organized shape and frontal density of fish schools. Ethology, 
114(3):245-254, 2008. 

M. O. Jackson, Social and Economic Networks. Princeton University Press, 2010. 

A. Jadbabaie, J. Lin, and A.S. Morse, Coordination of groups of mobile autonomous agents using nearest 
neighbor rules. IEEE Trans. Automat. Control, 48, 988-1001, 2003. 

T. Karper, A. Mellet and K. Trivisa, Hydrodynamic limit of the kinetic Cucker-Smale flocking model 
ArXiv.120b.6S31, 2012. 

U. Krause, A discrete nonlinear and non-autonomous model of consensus formation. Communications in 
difference equations, page 227-236, 2000. 

P. F. Lazarsfeld, R. K. Merton, Friendship as a Social Process: A Substantive and Methodological Analysis. 
In "Freedom and Control in Modern Society" (M. Berger, T. Abel, and C. H. Page, eds.) New York, Van 
Nostrand, 18-66, 1954. 

W. Li, Stability analysis of swarms with general topology IEEE Trans. Systems, Man, Cyber., Part B, 
38(4), 2008. 

H. Liu and E. Tadmor, Critical thresholds in convolution model for nonlinear conservation laws, SIAM J. 
Math. Anal. 33(4), 930-945, 2001. 

J.Lorenz, Continuous opinion dynamics of multidimensional allocation problems under bounded confidence. 
A survey. Internat. J. Modern Phys. C, 18, 1819-1838, 2007. 

M. McPherson, L. Smith-Lovin, and J. M Cook, Birds of a Feather: Homophily in Social Networks. Annual 
Review of Sociology, Vol. 27: 415-444 2001. 

N. Mecholsky, E. Ott and T. M. Antonsen, Obstacle and predator avoidance in a model for flocking Physica 

D. 239, 988-996, 2010. 

T. Mengesha and Q. Du, Analysis of a scalar peridynamic model with a sign changing kernel prteprint 

B. Mohar, The Laplacian spectrum of graphs, In " Graph Theory, Combinatorics, and Applications" ,(Y. 
Alavi, G. Chartrand, O. R. Oellermann and A. J. Schwenk, eds.), 2, John Wiley & Sons, 1991, pp. 871- 



HETEROPHILIOUS DYNAMICS ENHANCES CONSENSUS 



39 



[52] S. Motsch and E. Tadmor, A new model for self-organized dynamics and its flocking behavior. Journal of 

Statistical Physics, 144(5) :923-947, August 2011. 
[53] R. Olfati-Saber, J. A. Fax, and R.M. Murray, Consensus and cooperation in networked multi-agent systems. 

Proc. IEEE, 95, 215-233, 2007. 
[54] B. Piccoli and A. Tosin, Time-evolving measures and macroscopic modeling of pedestrian flow. Arch. Rat. 

Mech. Anal., 199, 707-738, 2011 
[55] C. W. Reynolds, Flocks, herds and schools: A distributed behavioral model. In ACM SIGGRAPH Computer 

Graphics, 21, 25-34, 1987. 

[56] S. Schochet and E. Tadmor, The regularized Chapman-Enskog expansion for scalar conservation laws, Arch. 

Rational Mech. Anal. 119, 95-107, 1992. 
[57] S. A. Silling, Reformulation of Elasticity Theory for Discontinuities and Long- Range Forces. J. Mech. 

Physics of Solids, 48, 175-209, 2000. 
[58] H. Spohn, Large Scale Dynamics of Interacting Particles. Texts and Monographs in Physics, Springer 1991. 
[59] E. Tadmor and C. Tan, Critical thresholds in compressible Eulerian dynamics with nonlocal viscosity 

preprint 

[60] J. Toner and Y. Tu, Flocks, herds, and schools. A quantitative theory of flocking. Physical Review E. 58, 
4828-4858 (1998). 

[61] C. M. Topaz and A. L. Bertozzi, Swarming patterns in a two-dimensional kinematic model for biological 

groups. SI AM J. Appl. Math. 65, 152-174 (2004). 
[62] G. Toscani, Kinetic models of opinion formation. Comm. Math. Sci., 4,481-496, 2006. 

[63] T. Vicsek, A. Czirok, E. Ben-Jacob, I. Cohen, and O. Shochet, Novel type of phase transition in a system 
of self-driven particles. Physical Review Letters, 75(6), 1226-1229, 1995. 

(Sebastien Motsch) 

Center for Scientific Computation And Mathematical Modeling (CSCAMM) 
University of Maryland, College Park, MD 20742 USA 

E-mail address: smotsch@cscamm.umd.edu 

URL: http://www.seb-motsch.com 

(Eitan Tadmor) 

Center for Scientific Computation And Mathematical Modeling (CSCAMM) 

AND 

Department of Mathematics, Institute for Physical Science and Technology 
University of Maryland, College Park, MD 20742 USA 

E-mail address: tadmor@cscamm.umd.edu 

URL: http : //www . cscamm . umd . edu/tadmor 



